**********************************************
*** PACs & January 6th ***********************
*** Kenneth M. Miller & Tanner Bates *********
**********************************************
* This code will reproduce all results using dta file:
* PACs and J6.dta

* Table 1 Models
* PAC %
probit anyreject biz_pct00 dw_one housevote trumpvote terms
estimates store m1
* PAC $ (in millions)
probit anyreject biz_mill dw_one housevote trumpvote terms
estimates store m2

estout m1 m2, cells(b(star fmt(3)) se(par fmt(3))) stats(r2 N, labels("Pseudo R2" "N"))

* Figure 1. Predicted probability with rugplot in R, see separate file:
* "Rug Plot for Jan 6.R"


* Supplemental Appendix

*A1. Business PAC totals
tabstat Total_Receipt, by(anyreject) stat(sum) format(%9.0f)
tabstat business, by(anyreject) stat(sum) format(%9.0f)


*A2. PAC percentages
tabstat business labor ideological other_pacs all_pacs, by(anyreject) stat(sum) format(%9.0f)


* A3. Additional models %:
* All PACs
probit anyreject pac_pct00 dw_one housevote trumpvote terms
estimates store m3
* Small donor percentage instead of PAC percentage
probit anyreject small_pct00 dw_one housevote trumpvote terms
estimates store m4
* Freedom Caucus
probit anyreject biz_pct00 freedom housevote trumpvote terms
estimates store m5
* Legislative effectiveness added to main model
probit anyreject biz_pct00 dw_one les housevote trumpvote terms
estimates store m6
* Roth rating instead of member percent
probit anyreject biz_pct00 dw_one safety trumpvote terms
estimates store m7
* Freshman indicator instead of terms in Congress
probit anyreject biz_pct00 dw_one housevote trumpvote freshman
estimates store m8
* Trump endorsement instead of terms in Congress
probit anyreject biz_pct00 dw_one housevote trumpvote endorse
estimates store m9

estout m3 m4 m5 m6 m7 m8 m9, cells(b(star fmt(3)) se(par fmt(3))) stats(r2 N, labels("Pseudo R2" "N"))


* A4. Additional models $:
* All PACs
probit anyreject pac_mill dw_one housevote trumpvote terms
estimates store m10
* Small donor percentage instead of PAC percentage
probit anyreject small_mill dw_one housevote trumpvote terms
estimates store m11
* Freedom Caucus
probit anyreject biz_mill freedom housevote trumpvote terms
estimates store m12
* Legislative effectiveness added to main model
probit anyreject biz_mill dw_one les housevote trumpvote terms
estimates store m13
* Roth rating instead of member percent
probit anyreject biz_mill dw_one safety trumpvote terms
estimates store m14
* Freshman indicator instead of terms in Congress
probit anyreject biz_mill dw_one housevote trumpvote freshman
estimates store m15
* Trump endorsement instead of terms in Congress
probit anyreject biz_mill dw_one housevote trumpvote endorse
estimates store m16

estout m10 m11 m12 m13 m14 m15 m16, cells(b(star fmt(3)) se(par fmt(3))) stats(r2 N, labels("Pseudo R\textsuperscript{2}" "N"))


* A5. Summary statistics
sum anyreject biz_pct00 biz_mill pac_pct00 pac_mill small_pct00 small_mill dw_one les freedom housevote safety trumpvote terms freshman endorse if anyreject==0 | anyreject==1

